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1.0  Introduction 


For  decades  the  U.S.  Navy  has  used  large  steel  barges  as  a  vital  part  of  their 
Military  Preposition  Force  (MPF)  offload  process.  The  MPF  is  a  fleet  of  cargo  ships 
prepositioned  at  strategic  locations  ready  for  deployment  when  called  upon.  For  any 
extended  military  action,  the  MPF  is  the  main  source  of  equipment  transport.  Many 
of  the  locations  to  which  the  U  S  military  may  be  called  into  action  do  not  have  deep 
water  facilities  that  can  accommodate  these  vessels.  As  a  result,  the  equipment  is 
offloaded  several  miles  off  shore  and  transported  to  the  beach  or  permanent/temporary 
piCT  via  barges  or  "causeways"  as  the  Navy  has  designated  them. 

The  causeways  currently  being  used  by  the  Navy  measure  21  foot  (ft)  wide,  90 
ft  long  and  have  a  5  ft  draft.  The  shape  is  very  close  to  a  rectangular  steel  box.  Due 
to  it’s  shallow  draft  and  particular  geometry,  the  causeways  quite  often  have  water  on 
the  deck  even  in  relatively  low  sea  states.  Hence  their  stability  is  sensitive  to  wave 
action.  To  improve  the  reliability  of  their  operation,  a  new  generation  of  causeways, 
in  which  the  dimensions  increase  to  25  ft  wide,  120  ft  long  and  a  draft  of  8  ft  is 
currently  under  conceptual  design  and  is  the  focus  of  this  analysis. 

1.1  Objective 

A  series  of  experiments  have  been  conducted  to  provide  data  to  calibrate  the 
prediction  capabilities  of  several  time,  frequency  and  probability  domain  computer 
programs  developed  by  the  Naval  Facilities  Engineering  Service  Center  (NFESC), 
Amphibious  Division.  Of  particular  interest  to  the  NFESC  engineers  is  a  study  of 
capsizing  when  the  causeways  are  subjected  to  high  sea  states  for  relatively  long 
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durations.  Simple  low  degree-of-freedom  dynamic  models  of  the  barge  motions  have 
been  developed  for  stochastic  analytical  and  numerical  studies  to  estimate  the 
probability  of  capsize.  The  motivation  behind  this  study  is  to  identify  suitable  system 
damping  parameters  for  the  single-degree-of-freedom  (SDOF)  model. 

1.2  Literature  Review 

The  motion  of  free-floating  vessels  in  open  seas  often  contain  complicated 
nonlinear  behaviors  due  to  large  displacement  and  hydrodynamic  effects.  In 
particular,  the  stability  of  roll  motion  of  barges  with  high  centers  of  gravity  is  of 
great  practical  importance.  Existing  stability  criteria  are  expressed  in  terms  of 
minimum  values  of  certain  key  features  of  the  righting  arm.  For  certain  classes  of 
vessels,  static  stability  standards  based  on  statistical  and  other  analyses  of  intact  static 
condition  are  sufficient  for  design  purpose,  and  can  give  a  qualitative  understanding  of 
the  stability  behavior  for  the  naval  architect  (Soliman  and  Thompson  1991).  Although 
this  curve  is  found  to  be  an  important  vessel  characteristic  in  assuring  safety 
(Falzarano  et  al.  1992),  other  vessel  properties  are  also  significant.  These  include 
hydrodynamic  and  viscous  roll  damping,  as  affected  by  the  wave  exciting  force. 

Because  the  cargo  is  carried  above  deck  on  a  causeway,  the  center  of  gravity 
of  the  entire  system  is  considerably  higher  in  proportion  to  that  of  a  conventional  ship 
which  carries  most  of  its  cargo  below  deck.  The  dynamics  of  a  high  center-of-gravity 
barge  rolling  in  seaways  can  become  difficult  to  model  and  predict  because  of  the 
highly  nonlinear  characteristics  encountered.  Accidental  capsize  of  log  barges  (which 


3 


have  very  similar  geometry  and  loading  conditions  as  those  considered  in  this  study) 
under  average  to  calm  sea  condition  has  been  well  documented  (McAllister  1995). 

1.3  Scope 

The  scope  of  this  study  is  to  model  the  barge  response  when  subjected  to  beam 
sea  conditions  as  a  SDOF  system  for  roll  motion.  The  measured  wave,  either  regular 
or  random,  is  input  into  a  computer  program  which  uses  a  fourth-order  Runge-Kutta 
integration  method  to  numerically  predict  roll  motion  response.  The  simulated  results 
are  compared  to  measured  data  to  determine  the  best  system  damping  parameters. 

Four  different  forms  of  the  damping  moment  of  the  SDOF  model  are  analyzed.  In 
each  form  all  coefficients  are  known  except  for  the  linear  and  nonlinear  damping 
parameters.  Each  form  uses  a  combination  of  relative  motion  "Morison"  damping 
and/or  linear  structural  damping.  Predicted  results  of  the  four  forms  of  the  damping 
momait  are  examined  to  determine  the  most  suitable  model.  A  sensitivity  study  on 
the  response  to  various  system  parameters  is  then  conducted  on  the  selected  form. 
Each  form  of  the  model  uses  a  thirteenth-order  polynomial  restoring  moment.  An 
analysis  is  also  conducted  to  determine  the  effects  of  using  lower  order  terms  to 
represent  the  restoring  moment.  Finally,  results  of  the  SDOF  model  from  this  report 
are  compared  with  a  3DOF  and  a  2DOF  model  subjected  to  the  same  wave 
environment  to  determine  which  model  more  accurately  predicts  roll  response. 
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2.0  System  Description 

The  system  considered  in  this  analysis  is  that  of  a  free  floating  barge  subjected 
to  regular  and  random  waves  in  open  seaways.  In  general,  a  six  degree-of-freedom 
(6DOF)  model  is  required  to  fully  specify  the  motions  of  the  barge.  Subjecting  a 
barge  to  beam  seas  would  be  considered  the  most  critical  scenario  for  barge  capsize. 
This  simplifies  the  equations  of  motion  because  under  ideal  beam  conditions,  the 
response  motions  for  a  barge  may  be  described  by  a  three  degree-of-freedom  (3DOF) 
model  with  heave,  sway  and  roll  represented  by 

wy  +  m^cos(|^)(y-<')  +  m^^sin(|^l)(y-i>)  +  Cjjfy-v) 

-  m^z  -  m(ZgCOs<f>)^ 

+  R^i<l>,z,ri,h)sm(h)  +  K^y  =  0 

mz  +  m  cos(^)(z-w)  +  sin(|4^|)(z-vi;)  + 

dy  \  ‘  (2.1) 

+  +  m4>y  “  m(zfos<f>)^  +  mg 

*  C0S(|J)  =  0 

dy  dy 

+  m(z  cos<j))^z  -  m(z  cos0)y  +  RJ<l>,z, v,^)cos(^) 

*  *  ^  dy  dy 

-  mgZgSin<j>  =  0 

(see  Bartel  1996  for  derivation  and  definitions).  Although  this  model  is  significantly 
simpler  than  a  large  number  of  time  domain  large  body  analysis  models  based  on 
potential  theory,  it  is  evident  that  these  equations  are  still  relatively  involved  and 
computationally  intensive  to  solve  numerically  compared  to  that  of  a  SDOF  model. 
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The  3D0F  model  based  on  these  equations  is  currently  being  analyzed  under  a 
separate  study  (Bartel  1996). 

2.1  Assumptions 

There  are  several  assumptions  made  in  the  development  and  application  of  the 
model  equations  for  this  analysis.  The  first  is  that  the  nominal  wave  length  is 
significantly  greater  than  the  beam  of  the  barge.  As  a  result,  the  wave  surface  can  be 
assumed  linear  across  the  beam.  Secondly,  wave  forces  and  moments  act  at  the 
center  of  gravity  in  lieu  of  integrating  pressure  over  the  surface.  Finally,  it  is 
assumed  that  roU  response  can  be  uncoupled  from  heave  and  sway,  therefore  the 
focus  will  be  only  on  the  pure  roll  motion  of  the  barge  modeled  by  a  SDOF  system. 

2.2  Equation  of  Motion 

Uncoupled  roll  motion  subjected  to  beam  sea  conditions,  assuming  long  wave 
excitation  and  deep  water  conditions,  is  modeled  by  the  equation 

Ij  +  Im4(^ -dv/dy)  +  C^^(^-dii/dy)  +  C^j,(^-drj/dy)  |  ^-dv/dy  \  (2.2) 

+  =  0 

where  744  is  the  rotational  inertia  (in  air)  of  the  vessel  about  the  roll  axis,  the 
hydrodynamic  added  inertia,  ^  the  roll  acceleration,  dfi/dy  the  wave  slope 
acceleration,  044^,  the  linear  damping  coefficient,  ^  the  roll  velocity,  drj/dy  the  wave 
slope  velocity,  €44^  the  nonlinear  damping  coefficient,  and  R44(<i>,i],dr]/dy)  the 
nonlinear  restoring  moment. 
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2.3  Description  of  Wave  Field 

Based  on  linear  wave  theory,  the  velocity  potential  and  related  wave  profile, 
time  derivatives  and  partial  spatial  derivatives  are  as  follows: 


$ 

V 

V 

V 


=  -  <«) 
(j)  cosh(^) 

"  "  ^cos(^y  -  wr) 

g  ot 

=  ^  =  w^sin(^y  -  wf) 

=  ^  =  -aPAcos(/^-wt) 


(2.3) 


1}'  =  ^  =-  ^^sin(^y  -  at) 
dy 

ri'  =  ^  =  akAcos(ky  -  at) 

^'  =  —  =  ct)^A:/4sin(^y  -  at) 
dy 


If  deep  water  wave  conditions  are  applied  and  water  particle  characteristics  are 


considered  at  the  surface,  the  wave  number,  k,  can  be  represented  by 


k= 


a‘- 


g 


(2.4) 
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and  the  resulting  wave  field  characteristics  reduce  to 

ri  =  ^cos(^  -  o)t) 
i]  =  («i^4sin(^  -  (at) 
ri  = 


2.4  Inertia  and  Added  Inertia  Coefficients 

For  a  given  barge,  the  rotational  inertia,  can  be  analytically  computed. 

The  hydrodynamic  added  inertia,  1^,  depends  on  the  frequency  of  the  barge  roll 
motion.  The  value  of  the  added  inertia  under  the  range  of  frequencies  of  the  roll 
motion  response  observed  in  this  study  is  relatively  constant.  Therefore,  a  constant 
added  inertia  coefficient  is  adopted  for  convenience  of  analysis.  The  first  two  terms 
of  equation  2.2  represent  the  inertial  moment  of  the  barge  in  air,  lji>,  and  the  relative 
motion  added  inertia  moment,  lA^lQ^-^l^y),  respectively.  The  term  (^dfi/dy)  is  the 
relative  motion  between  the  roll  acceleration  of  the  barge  and  the  wave  slope 
acceleration  resulting  from  wave  excitation.  For  this  analysis,  the  value  of  144  = 
2.161046e6  slug-ft^  (Paulling  1995).  With  this  information  and  knowing  the  natural 
period  of  the  barge  is  approximately  5  seconds,  the  added  inertia  is  calculated  to  be 
Ia44  =  1.3725e6  slug-f^. 
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2.5  Damping  Moment 

In  equation  2.2,  the  damping  moment  resulting  from  radiation  of  waves  from 
the  barge  (corresponding  to  linear  energy  dissipation  in  potential  theory)  is 
represented  by  the  linear  term,  C44iX^-dii/dy).  The  effects  of  damping  resulting  from 
flow  separation  and  turbulence  are  represented  by  the  nonlinear  Morison  term, 

Q4n(^  ■  dr}/dy)  \  ^  -  dri/dy  \ .  €44^,  and  €44^,  are  the  linear  and  nonlinear  damping 
coefficients  expressed  in  the  following  convenient  form: 


2  * 


'UL 


C  = 

'-'44Ar 


2  * 


(2.6) 


where  «„  is  the  "natural  frequency"  of  the  barge  (determined  by  small  amplitude  free 
vibration  tests)  and  is  equal  to  approximately  5  seconds  (Yim  et  al.  1995).  and 
are  the  damping  parameters  associated  with  the  linear  and  nonlinear  damping  terms 
(to  be  determined  by  system  identification  and/or  numerical  iteration).  The  terms 
have  been  normalized  by  the  sum  of  the  inertia  and  added  inertia  terms,  Ij  =1^4  + 

7^.  The  term  (0  -  drildy)  is  the  relative  motion  between  the  roll  velocity  of  the  barge 
and  the  wave  slope  velocity  resulting  from  wave  excitation. 


2.6  Roll  R^oring  Moment 

The  restoring  moment  term,  R44i<l>,ri,dr]fdy),  represents  the  moment  created 
from  buoyancy  and  gravity  forces  acting  on  the  barge  under  static  conditions.  The 
restoring  moment  is  a  thirteenth-order  (13th-order)  polynomial  made  up  of 
coefficients  of  the  righting  moment  curve  multiplied  by  the  corresponding  power  of 
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the  relative  motion  between  the  roll  displacement  and  the  wave  slope  as  shown  in  the 
following  equation: 


B,,(<l>-dv/dyy 

4 


^  B,^i<l>-dvldyf 


+ 


^11,5 

B,,{<t>-dvldyy\ 

It 


B,^{<i>-drildy)\ 

h 

B,,(<{>-dnldyr 

4 


(2.7) 


The  coefficients  Bi  5,  B3  5,  etc.  are  established  by  fitting  a  thirteenth-order  polynomial 
to  the  lighting-moment  curve  shown  in  Figure  2.1  (Yim  et  al.  1995).  The  term  is 
normalized  by  the  total  inertia,  Ij-.  Although  a  lower-order  polynomial  can  be  used 
with  sufficient  accuracy  for  small  angle  motions,  the  additional  time  and  effort 
associated  with  using  a  thirteenth-order  polynomial  is  minimal,  and  in  the  case  of 
large  angle  motions,  the  higher  order  terms  are  essential. 


2.7  Various  Forms  of  Roll  Equation  of  Motion 

Several  combinations  of  the  terms  representing  the  damping  moment  in  the 
equation  of  motion  are  examined  here.  In  addition  to  equation  2.2  which  uses  relative 
motions  in  both  the  linear  and  nonlinear  damping  terms,  models  which  approximately 
account  for  and  neglect  different  physical  effects  wiU  be  analyzed.  One  combination 
assumes  the  magnitude  of  the  wave-slope  velocity  is  negligible  compared  to  the  roll 
velocity  (e.g.  near  roll  resonance).  As  a  result  only  the  roll  velocity  component  will 
be  considered  in  the  linear  and  nonlinear  damping  terms.  This  is  modeled  by 
equation  2.8: 


Righting  Moment  (R-lb) 


Ron  (deg) 


Figure  2.1  13th  Order  Roll  Restoring  Moment  Curve.  -  numerical  (Paulling  1995) 

--  polynomial  fit  (Yim  et  al.  1995).  B,3_s  =  5.605 le6,  B„5  =-1.1657e6, 

B9,5=  1.7258e6,  9.6915e6,  B5,5=  -9.3298e6,  B3,5=  3.3999e6,  B,  5=  -4.4569e6 
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+  R^(^,ri,drildy)  =  0 

Equation  2.9  below  assumes  that  the  linear  damping  coefficient  is  dominated 
by  the  barge  characteristics  (damping  induced  by  the  presence  of  mooring  lines  in  the 
e7q)eriments  which  is  not  modeled  explicitly  in  this  study).  The  relative  motion 
between  the  barge  and  wave  slope  velocity  is  accounted  for  only  in  the  nonlinear 
damping  term. 


+  C44^<^  +  -dr]/dy)  \  ^  -  dij/dy  \ 

+  R^{<t>,i)Mldy)  =  0 

The  final  form  of  the  model  to  be  analyzed  is  equation  2.10: 


(2.9) 


4m (^ “ 3^/3y)  +  C44^(<^ - biildy)  \^-diilBy\ 
+  R^^(<l>,ri,dr}ldy)  =  0 


(2.10) 


This  model  assumes  that  damping  effects  are  strictly  nonlinear  and  relative  motion 
between  the  wave  slope  velocity  and  barge  roll  velocity  significantly  influences  the 
roll  response. 
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3.0  Model  Tests  and  Identification  of  Damping  Parameters 

The  four  different  forms  of  the  model  for  roll  motions  are  analyzed  to 
determine  which  form  best  predicts  roll  response  when  the  barge  is  subjected  to  beam 
sea  conditions.  In  each  case,  the  prototype  measured  wave  is  sampled  at  2  Hertz 
(Hz).  The  data  is  then  filtered  so  that  frequencies  above  0.25  Hz  are  disregarded 
because  they  do  not  contribute  to  the  response  motion.  Wave  profile  and  water 
particle  velocity  as  well  as  wave  slope,  wave  slope  velocity  and  wave  slope 
acceleration  are  then  numerically  computed.  Regular  wave  cases  analyzed  are  that  of 
the  5,  6,  8,  and  10  sec  waves  with  a  desired  6  ft  wave  height.  In  addition,  two 
random  cases  are  modeled;  a  Bretschneider  spectrum  of  6  ft  significant  wave  height 
and  8  sec  dominant  wave  period,  and  a  white  noise  spectrum  of  6  ft  significant  wave 
height  and  wave  periods  ranging  from  4  to  20  secs. 

In  each  of  the  forms,  all  terms  are  known  except  for  the  linear  and  nonlinear 
damping  parameters,  1^,  and  used  in  the  damping  coefficients,  €44^,  and  €44;^.  An 
initial  estimate  for  the  damping  parameters  is  obtained  through  system  identification 
procedures.  By  comparing  the  standard  deviation  of  predicted  roll  response  to 
measured  roll  response,  an  iterative  approach  is  then  used  to  identify  the  most 
appropriate  damping  parameters  for  the  model.  The  inertia  coefficients  and  the 
nonlinear  restoring  moment  coefficients,  constant  for  all  the  wave  cases,  are 
determined  through  experimental  measurements  or  computer  analysis.  The  other 
coefficients  are  determined  through  the  measured  data  for  each  particular  wave 


excitation  test. 
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3.1  Model  Test  Procedure 

Regular  and  random  wave  model  tests  have  been  conducted  on  a  single  barge 
to  simulate  response  in  opai  seas.  The  facility  used  to  conduct  the  tests  was  a  48  ft 
wide,  250  ft  long  basin  which  was  cjq)able  of  generating  long  crested  waves  up  to  3  ft 
in  height.  A  1/16  scale  model  of  the  120  ft  long,  25  ft  wide  and  8  ft  deep  full  scale 
barge  was  used.  For  this  analysis,  only  results  of  the  barge  subjected  to  beam  sea 
conditions  are  examined.  Although  the  objective  of  the  tests  is  to  simulate  motions  of 
the  barge  in  open  seas,  mooring  lines  w&re  attached  to  the  model  to  prevent  it  from 
drifting  down  the  basin.  The  Iragth,  weight  and  sag  of  the  mooring  lines  are 
designed  to  ensure  the  natural  period  of  the  mooring  system  would  be  much  larger 
than  the  wave  excitation  pmod. 

The  particular  wave  conditions  generated  for  this  report  were  4  regular  wave 
and  2  random  wave  tests.  Regular  waves  of  approximately  6  ft  wave  height  (full 
scale)  were  conducted  for  wave  periods  of  5,  6,  8  and  10  secs.  The  random  wave 
tests  were  for  a  Bretschneider  spectrum  of  6  ft  significant  wave  height  with  and  8  sec 
dominant  wave  period,  and  for  a  band  limited  "white  noise"  test  of  6  ft  significant 
wave  height  with  wave  periods  ranging  from  4  to  20  secs.  Model  tests  were  run  until 
steady  state  had  been  reached  before  data  collection  began.  Only  1024  secs  (17 
minutes)  of  data  was  collected  for  each  run  in  order  to  avoid  significant  energy 
accumulation  from  wave  reflection,  re-reflection  and  other  "tank  noise"  which  might 
build  up  over  time  and  could  corrupt  the  wave  excitation. 
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3.2  Description  of  Analytical  Prediction  Procedure 

Test  SB26,  the  6  ft,  5  sec  regular  wave,  will  be  used  as  an  example  of  the 
procedure  followed  to  establish  damping  coefficients  and  simulated  roll  response. 

The  entire  time  series  of  the  measured  wave  input  is  examined  to  confirm  the  desired 
wave  height  and  period  was  consistent  with  the  target  spectrum.  Figure  3.  la  shows 
the  1024  sec  time  series  used  in  this  analysis.  The  measured  data  is  then  subjected  to 
a  low  pass  filter  with  cutoff  frequency  of  0.25  Hz  after  digital  sampling.  The  filtered 
wave  is  then  displayed  as  shown  in  Figure  3.1b,  this  is  the  wave  data  input  into  the 
simulation  program.  To  ensure  the  input  wave  for  simulation  is  correct,  the  spectral 
densities  of  the  measured  and  filtered  wave  are  compared  as  shown  in  Figure  3.2. 

The  damping  parameters  for  the  linear  and  nonlinear  damping  coefficients  are  input 
into  the  model  equation  and  the  simulation  is  executed.  The  results  of  the  measured 
roll  response.  Figure  3.1c,  are  then  compared  to  those  of  the  simulated  roll  response. 
Figure  3.  Id.  A  small  section  of  the  time  series  is  compared  as  well,  as  shown  in 
Figure  3.3.  The  standard  deviation  of  the  measured  and  simulated  data  is  then 
compared  to  ensure  the  most  accurate  match  of  roll  response.  In  the  random  wave 
tests,  the  above  procedure  is  performed  and  an  additional  comparison  of  the  spectral 
and  probability  densities  of  roll  response  are  made  to  ensure  the  most  accurate  match 
of  roll  response. 


roll  (degrees)  roll  (degrees) 
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(a)  Measured  Wave  Profile 


(b)  Filtered  Wave  Profile 


(c)  Measured  Roll 


(d)  Predicted  Roll  time  (seconds) 


Figure  3.1  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Regular  Wave  H=  6.6  Ft,  T=  5  Sec  (SB26) 


Gxx 
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roll  (degrees)  and  wave  profile  (feet) 
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3.3  Form  1:  Relative  Motion  Damping  Moment 

The  first  model  to  be  analyzed  is  that  of  equation  2.2: 

-drildy)  +  C44j.(0  -diildy)  +  C^yi^-drildy)  |  ^  -di)/dy  [  «  2\ 

*  R^(4>,v,dvfdy)  =0 

Recall  this  form  assumes  both  the  linear  (radiation  damping)  and  nonlinear  (Morison 
type  local  flow  separation)  have  significant  influence  on  the  roll  motion  response. 

The  resulting  linear  and  nonlinear  damping  coefficients  for  this  model  are  listed  in 
Table  3.1: 


Wave  Test 

Linear  Damping 
Paramet^ 

Nonlinear  Damping 
Parameter 

H=6.6  ft,  T=5  sec 
(SB26) 

0.19 

0.19 

H=6  ft,  T=6  sec 
(SB27) 

0.11 

0.11 

H=7.2  ft,  T=8  sec 
(SB29) 

0.80 

0.20 

H=5.7  ft,  T=10  sec 
(SB30) 

0.03 

0.07 

Bietschneider 

Spectrum  (SB25) 

0.04 

0.04 

White  noise 

Spectrum  (SB33) 

0.08 

0.06 

Table  3.1  linear  and  Nonlinear  Damping  Parameters  for  Form  1,  Equation  2.2. 
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3.3.1  Regular  Wave,  H=  6  Ft,  T=  5  Sec  (SB26) 

Figure  3.4  represents  the  roll  response  of  the  model  barge  subjected  to  a  6.6 
ft  wave  with  a  5  sec  period.  The  maximum  roll  response  is  approximately  13 
degrees.  The  model  indicates  a  linear  damping  parameter  of  0.19  and  a  nonlinear 
danqnng  parameter  of  0.19.  The  amplitude  of  the  measured  and  predicted  re^nse  is 
nearly  identical.  However,  the  predicted  roll  re^nse  is  out  of  phase  by 
approximately  1  sec  (1/5  period)  with  the  measured  data. 

3.3  J:  R^ular  Wave,  H=  6  Ft,  Ts=  6  Sec  (SB27) 

Figure  3.5  rqiresents  the  results  of  the  model  barge  subjected  to  a  6  ft  regular 
wave  widi  a  6  sec  wave  period.  The  maximum  roll  response  is  11.5  degrees.  The 
model  indicates  both  the  linear  and  nonlinear  damping  parameters  are  0.11.  Figure 
3.5  shows  that  amplitude  of  the  predicted  roll  response  matches  the  measured  data 
well.  As  in  the  5  sec  wave  case,  the  predicted  roll  response  is  approximately  1  sec 
(1/6  period)  out  of  phase  with  the  measured  data. 

3.3.3  R^ular  Wave,  H-  7.2  Ft,  T=  8  Sec  (SB29) 

Figure  3.6  represents  the  results  of  the  model  barge  subjected  to  a  7.2  ft 
regular  wave  with  an  8  s^  wave  period.  The  damping  parameters  established  for  this 
wave  test  are  0.80  for  linear  damping  and  0.20  for  nonlinear  damping.  Note  that  the 
linear  damping  parameter  is  agnificandy  higher  than  the  two  previous  test  cases.  The 
maximum  roll  re^nse  of  the  measured  data  is  4.1  degrees  and  is  well  matched  by 


time  (seconds) 

Figure  3.4  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Regular  Wave,  H=  6.6  Ft,  T=  5  Sec  (SB26,  eqn  2.2) 

-  measured  wave,  -.  predicted  wave,  ~  measured  roll,  ..  predicted  roll. 


time  (seconds) 

Figure  3.5  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Regular  Wave,  H=  6  Ft,  T=  6  Sec  (SB27,  eqn  2.2) 

-  measured  wave,  -.  predicted  wave,  —  measured  roll,  ..  predicted  roll. 


time  (seconds) 

Figure  3.6  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Regular  Wave  H=  7.2  Ft,  T=  8  Sec  (SB29,  eqn  2.2) 

-  measured  wave,  -.  predicted  wave,  -  measured  roll,  ..  predicted  roll. 
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the  predicted  reqwnse  of  the  model.  This  form  of  the  model  is  not  able  to 
consistaitly  match  the  phase  of  the  response.  The  predicted  data  appears  to  be 
approximately  1  sec  (1/8  period)  out  of  phase  (at  the  peak)  with  the  measured  data. 

3.3.4  Regular  Wave,  H=  5.7  Ft,  T=  10  Sec  (SB30) 

Test  SB30,  in  which  the  excitation  is  a  5.7  ft  wave  with  a  10  sec  period,  proves 
to  be  the  most  difficult  reqwnse  to  match  by  any  form  of  the  model.  Because  the 
natural  period  of  the  barge  is  approximately  5  secs,  which  is  a  multiple  of  the  wave 
exdtation  period,  the  response  contains  a  siq)erharmonic. 

Figure  3.7  represents  the  best  comparison  of  standard  deviation  and 
superharmonic  response  which  could  be  obtained  with  the  givai  form  of  the  model 
equation  (equation  2.2).  If  the  damping  parameter  is  increased,  the  superharmonic 
characteristics  are  eliminated.  If  the  damping  parameter  is  decreased,  the  gimniatpri 
response  amplitude  is  much  greater  than  the  measured  reqwnse  amplitude.  The 
damping  parameters  used  in  this  figure  are  0.03  for  the  linear  coefficient  and  0.07  for 
the  nonlinear  coefficirat.  The  standard  deviation  for  the  measured  re^nse  is 
calculated  to  be  1.84.  The  standard  deviation  of  the  simulated  response  is  1.90. 

Figure  3.8  shows  a  good  comparison  of  the  measured  and  simulated  spectral  dengities 
for  roll  re^nse. 


0  200  400  600  800  1000 

time  (seconds) 

Figure  3.9  1024  Sec  Time  Series  of  Predicted  Roll,  H=  5.7  Ft,  T=  10  Sec 
(SB30,  eqn  2.2) 
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Figure  3.9  provides  a  rq)reseiitation  of  predicted  roll  response  for  the  entire  1024  sec 
time  s^es.  It  is  evident  that  the  supeiiiarmonic  does  not  occur  in  all  locations  with 
the  same  amplitude.  A  parametric  study  presented  in  a  later  section  of  this  report  will 
provide  a  better  understanding  of  the  effects  of  varying  the  linear  and  nonlinear 
damping  parameters. 

3.3.5  Bretsdineider  Spectrum,  Hs=  6  Ft,  Tp=  8  Sec  (SB25) 

Figure  3.10a  -  3.10d  rq)resent  the  measured  and  filtered  wave  input  and  the 
resulting  measured  roll  response  and  predicted  roll  response  (to  filtered  wave 
recitation)  fOT  a  Bretschneider  wave  spectrum.  As  shown  in  Figure  3.11,  the  spectral 
denaties  of  die  measured  and  filtered  wave  excitation  input  are  compared  to  ensure 
the  proper  input  is  being  used  for  prediction.  The  damping  parameters  used  for 
prediction  are  equal  to  0.04  for  both  the  linear  and  the  nonlinear  terms.  A  visual 
comparison  of  Figures  3.10c  and  3.10d  indicates  the  model  provides  an  accurate 
prediction.  The  standard  deviation  for  the  predicted  roll  response  is  3.95  which 
compared  favorably  to  the  standard  deviation  of  3.95  for  the  measured  data.  Figure 
3.12  indicates  the  qiectral  doisities  are  a  very  close  match.  The  final  comparison 
made  between  measured  and  simulated  data  is  a  histogram  of  roll  reqxmse.  Figures 
3.13a  and  b  confirm  that  the  model  provides  an  accurate  representation  of  roU 
response  motion. 
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Figure  3.10  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Bretschneider  Spectrum,  Hs=  6  Ft,  Tp=  8  Sec  (SB25,  eqn  2.2) 
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3.3.6  White  Noise  Excitation,  Hs=  6  Ft  (SB33) 

Figures  3.14a  through  3.14d  r^resent  the  measured  and  predicted  wave  input 
and  measured  and  predicted  roll  response  for  a  white  noise  excitation  with  wave 
periods  of  4  to  20  secs  and  a  6  ft  significant  wave  height.  A  comparison  of  the 
spectral  densities,  Figure  3.15,  confirms  that  the  simulated  wave  input  m^itche-s  the 
measured  wave  excitation.  The  linear  and  nonlinear  damping  parameters  used  in  the 
mcxlel  are  0.08  and  0.06  respectively.  A  visual  comparison  of  Figures  3.14c  and 
3.14d  indicate  the  model  provides  a  good  simulation  of  roll  response.  The  standard 
deviation  of  predicted  roll  response,  5.30,  compared  well  with  the  measured  roU 
re^nse  which  is  5.29.  A  comparison  of  spectral  densities.  Figures  3.16,  and 
histograms.  Figures  3.17a  and  3.17b,  confirm  that  measured  and  predicted  roll 
respcmse  compare  quite  well.  Note  that  tiie  there  is  a  significant  increase  in 
probability  mass  at  the  extreme  values.  This  behavior  is  typical  of  nonlinear  systems 
with  a  softening  stifihiess  (which  is  the  case  for  barge  roll  motion). 

3.4  Fonn  2:  Relative  Motion  N^lected  in  Damping  Moment 

The  second  form  of  damping  representation  to  be  examined  is  the  following 

^44^  +  7^(0-^)  +  C44J.0  +  ^2.8) 

dy 

In  this  form,  the  relative  motion  between  the  wave  slope  velocity  and  the  barge  roll 
velocity  is  n^lected.  As  a  result,  only  the  barge  roll  velocity  contributes  to  the 
equation.  The  following  results  will  show  this  model  displays  an  equal  ability  to 
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predict  re^nse  amplitude  and  a  marked  improvement  in  matching  the  phase  of  the 
measured  roll  response.  Table  3.2  lists  the  resulting  damping  parameters  for  this 
model. 


Test  Case 

Linear  Damping 
Paramet^ 

Nonlinear  Damping 
Parameter 

H=6.6  ft,  T=5  sec 
(SB26) 

0.09 

0.07 

H=6  ft,  T=6  sec 
(SB27) 

0.10 

0.07 

H=7.2  ft,  T=8  sec 
(SB29) 

0.33 

0.40 

H=5.7  ft,  T=10  sec 
(SB30) 

0.03 

0.07 

Bretschneider  Spectrum 
(SB25) 

0.04 

0.04 

White  Noise  Spectrum 
(SB33) 

0.07 

0.08 

Table  3.2  linear  and  Nonlinear  Damping  Parameters  for  Form  2,  Equation  2.8 

3.4.1  Regular  Wave,  H=s  6.6  Ft,  T=  5  Sec  (SB26) 

Figure  3.18  displays  die  results  of  the  model  barge  subjected  to  a  6.6  ft  wave 
with  a  5  sec  period  considering  only  barge  roll  velocity.  The  damping  parameters 
used  to  achieve  these  results  are  0.09  for  the  linear  term  and  0.07  for  the  nonlinear 
term.  The  predicted  amplitude  closely  matches  the  measured  data.  In  addition,  the 
predicted  roU  response  motion  is  nearly  identical  in  phase  to  that  of  the  measured 
data. 
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Figure  3.18  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
RoU,  Regular  Wave,  H=  6.6  Ft,  T=  5  Sec  (SB26,  eqn  2.8) 

-  measured  wave,  -.  predicted  wave,  —  measured  roll,  ..  predicted  roll. 
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3.4.2  Regular  Wave,  H=  6  R,  T=  6  Sec  (SB27) 

Figure  3.19  r^esents  the  barge  response  when  subjected  to  a  6  ft,  6  sec  wave 
considering  only  barge  roll  velocity  in  the  damping  term.  The  damping  parameters 
used  to  achieve  these  results  are  0.10  for  the  linear  term  and  0.07  for  the  nonlinear 
term.  The  response  amplitude  adequately  matches  the  measured  reqx>nse  of  11.5 
degrees.  However,  the  predicted  roll  response  is  approximately  0.5  sec  (1/12  period) 
out  of  phase  with  the  measured  data.  Although  not  an  exact  match,  this  is  an 
improvement  over  the  phase  difference  of  approximately  1  sec  (1/6  period)  produced 
in  the  relative  motion  damping  form. 

3.4.3  R^ular  Wave,  H=  7.2  R,  H=  8  Sec  (SB29) 

Figure  3.20  rq)resents  the  barge  reqwnse  when  subjected  to  a  7.2  ft,  8  sec 
period  wave  considering  only  the  barge  roll  velodty  in  the  damping  term  for 
predicted  r^ponse.  The  damping  parameters  used  to  match  the  measured  roU 
response  of  4.1  degrees  are  0.33  for  the  linear  coeffidait  and  0.40  for  the  nonlinear 
coefficient.  As  with  the  relative  motion  damping  case,  these  coefficients  are 
significantly  higher  than  the  previous  test  cases.  Furtho-  study  should  be  conducted  to 
evaluate  this  phenomoion.  This  form  of  the  model  is  an  improvement  over  the 
relative  motion  damping  form  because  it  more  adequately  matches  the  measured 
behavior,  particularly  in  the  phase  (at  the  peaks). 
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Figure  3.19  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
RoU,  Regular  Wave,  H=  6  Ft,  T=  6  Sec  (SB27,  eqn  2.8) 

-  measured  wave,  predicted  wave,  --  measured  roll,  ..  predicted  roll. 
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Figure  3.20  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
RoU,  Regular  Wave,  H=  7.2  Ft,  T=  8  Sec  (SB29,  eqn  2.8) 

-  measured  wave,  -.  predicted  wave,  —  measured  roll,  ..  predicted  roll. 
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3.4.4  Regular  Wave,  H=  5.7  Ft,  T=  10  Sec  (SB30) 

Figure  3.21  rq)resents  the  roll  response  of  the  barge  subjected  to  5.7  ft  regular 
waves  with  a  10  sec  period,  considering  only  barge  roll  velocity  in  the  damping  term 
of  the  predicted  re^nse.  As  with  the  relative  motion  case,  it  is  difficult  to  match 
response  amplitude  with  the  siq)erharmonic  effect.  A  damping  parameter  of  0.03  for 
the  linear  damping  term  and  0.07  for  the  nonlinear  term  provides  the  closest  match  to 
the  measured  roll  response  of  2  degrees.  The  standard  deviation  of  the  measured  data 
is  1.84  compared  with  1.87  for  the  simulated  roll  response.  This  is  a  small 
improvemoit  over  the  prediction  using  relative  motion  damping  which  resulted  in  a 
standard  deviation  of  1.88.  Figure  3.22  shows  that  the  spectral  densities  of  the 
measured  and  simulated  roll  response  compare  favorably  for  this  form  of  the  model. 

3.4.5  Bretsdineider  Spectrum,  Hs=  6  Ft,  Tp=  8  Sec  (SB25) 

Figures  3.23a  -  3.23d  di^lay  the  measured  and  simulated  wave  exdtation  and 
measured  and  predicted  roll  req>onse.  The  damping  parameters  used  in  the  model  are 
0.04  fw  both  the  linear  and  nonlinear  damping  terms.  Figures  3.23c  and  3.23d 
indicate  a  good  visual  comparison  between  the  measured  and  simulated  reqwnses. 

The  standard  deviation  of  3.96  for  measured  data  compared  well  to  a  predicted 
standard  deviation  of  3.96.  The  energy  under  measured  and  simidated  spectral 
density  curves  are  a  close  match  although  there  is  a  slight  shift  in  peak  frequencies,  as 
shown  in  Figure  3.24.  A  final  comparison  of  probability  densities.  Figures  3.25a  and 
3.25b,  confirm  that  this  form  of  the  model  provides  a  good  simulation  for  the  random 
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Figure  3.21  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  Regular  Wave,  H=  5.7  Ft,  T=  10  Sec  (SB30,  eqn  2.8) 

-  measured  wave,  -.  predicted  wave,  —  measured  roll,  ..  pre^cted  roll. 
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3.5.6  IfVhite  Noise  Excitation,  Hs=  6  Ft  (SB33) 

Figures  3.26a  -  3.26d  diq)lay  the  measured  and  filtered  wave  and  predicted  roll 
response  for  a  white  noise  wave  exdtation  with  wave  periods  of  4  to  20  secs  and  a 
wave  hdght  of  6  ft.  The  damping  parametors  used  in  this  model  are  0.07  for  the 
linear  term  and  0.08  for  the  nonlinear  term.  The  resulting  standard  deviation  for 
predicted  roll  response  is  5.30  which  is  very  close  to  the  measured  req>onse  standard 
deviation  of  5.29.  A  comparison  of  q)ectral  densities,  Figure  3.27,  and  probability 
draisities.  Figures  3.28a  and  3.28b,  indicate  the  model  provides  an  accurate  prediction 
of  roll  response.  Again  a  significant  build-up  in  the  extreme  values  is  obsmved. 

3.5  Alternative  Fonns  Considered 

Other  forms  of  die  damping  moment  studied  are  those  represented  by  equations 
2.9  and  2.10: 

+  C^(i-dii/dy)  \  ^-drifdy  \  n  10) 

+  R^(4>,ri,dv/dy)  =0 

These  forms  are  capable  of  simulating  the  response  amplitude,  however  they  are  not 
as  accurate  in  duplicating  the  measured  phase.  Figure  3.29a  -  3.29d  is  a  comparison 
of  all  four  forms  for  test  SB26.  The  two  forms  which  use  relative  motion  in  the 
damping  term  have  the  greatest  phase  error.  This  indicates  the  system  is  inertia 
dominated  and  that  linear  (structural)  damping  makes  a  significant  contribution  to  the 
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Figure  3.26  Comparison  of  Measured  and  Filtered  Wave,  Measured  and  Predicted 
Roll,  White  Noise  Spectrum,  Hs=  6  Ft  (SB33,  eqn  2.8) 
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(a)  Measured  Roll  Amplitude  (degrees) 


(b)  Predicted  Roll  Amplitude  (degrees) 

Figure  3.28  Histograms  of  Measured  and  Predicted  Roll,  White  Noise  Spectrum 
(SB33,  eqn  2.2) 
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(a)  Form  1,  Equation  2.2 


(d)  Form  4,  Equation  2.10  time  (seconds) 

Figure  3.29  Comparison  of  Roll  Responses  Using  Four  Different  Forms  (SB26) 

-  measured,  --  predicted 


52 


re^nse.  Hie  form  which  uses  linear  damping  with  roll  velocity  only  and  relative 
motion  Morison  nonlinear  damping  provides  an  improvement  in  the  phase  error,  but 
is  still  slightly  out  of  phase.  Figure  3.29b  represents  the  model  which  uses  strictly 
structural  damping.  It  is  clear  this  model  provides  the  best  match.  This  trend  is 
evident  for  the  6  and  8  sec  wave  cases  as  well  as  shown  in  Figures  3.30  and  3.31. 

The  10  sec  wave  test  case  displayed  the  same  phase  shifted  results  for  all  model  forms 
as  shown  in  Figure  3.32. 

In  the  random  wave  test  cases,  each  model  is  equally  successful  at  simulating 
the  measured  refuse.  Equation  2.10,  which  uses  strictly  Morison  damping, 
required  a  significant  increase  in  the  damping  parameter  for  both  random  wave  cases. 
For  example  with  test  SB25,  the  Bretschneider  wave  exdtation,  the  damping 
parameters  used  for  equation  2.2,  2.8,  and  2.9  are  very  close  to  0.04  for  botii  the 
linear  and  nonlinear  damping  terms.  The  damping  parameter  for  equation  2.10, 
which  used  only  Morison  damping,  is  0.42.  Similar  results  are  obtained  with  test 
SB33,  the  white  noise  exdtation.  For  equations  2.2,  2.8  and  2.9,  the  damping 
parameters  are  in  the  range  of  0.07  to  0.08,  whCTe  as  the  damping  parameter  for 
model  equation  2.10  is  0.56.  Tables  3.3  and  3.4  summarize  the  resulting  damping 
parameters  for  equations  2.9  and  2.10. 

Form  4  of  the  SDOF  model  considered  only  relative  motion  Morison  damping. 
It  should  be  noted  that  all  the  damping  parameters  listed  in  Table  3.4  are  extremely 
high  for  the  regular  wave  cases.  These  ratios  are  unrealistic  and  indicate  that  Form  4 
is  not  an  acceptable  form  of  the  model. 
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(a)  Form  1,  Equation  2.2 


(d)  Form  4,  Equation  2.10  time  (seconds) 

Figure  3.31  Comparison  of  Roll  Responses  Using  Four  Different  Forms  (SB29) 
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(a)  Fonn  1,  Equation  2.2 


(b)  Form  2,  Equation  2.8 
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(c)  Form  3,  Equation  2.9 


(d)  Form  4,  Equation  2.10  time  (seconds) 

Figure  3.32  Comparison  of  Roll  Responses  Using  Four  Different  Forms  (SB30) 
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Wave  Test 

Linear  Damping 
Parameter 

Nonlinear  Damping 
Parameter 

H=6.6  ft,  T=  5  sec 
(SB26) 

0.09 

0.14 

H=6  ft,  T=6  sec 
(SB27) 

0.10 

0.11 

H=7.2  ft,  T=8  sec 
(SB29) 

0.34 

0.34 

H=5.7  ft,  T=10  sec 
(SB30) 

0.03 

0.07 

Bretschneider  Spectrum 
(SB25) 

0.04 

0.04 

White  Noise  Spectrum 
(SB33) 

0.07 

0.08 

Table  3.3  Linear  and  Nonlinear  Damping  Parameters  for  Form  3,  Equation  2.9 


Wave  Test 

Nonlinear  Damping 
Parameter 

H=6.6  ft,  T=5  sec 
(SB26) 

1.07 

H=6  ft,  T=6  sec 
(SB27) 

1.25 

H=7.2  ft,  T=8  sec 
(SB29) 

88.00 

H=5.7  ft,  T=10sec 
(SB30) 

1.25 

Bretschneider  Spectrum 
(SB25) 

0.42 

White  Noise  Spectrum 
(SB33) 

0.56 

Table  3.4  Linear  and  Nonlinear  Damping  Parameters  for  Form  4,  Equation  2.10 
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4.0  Sensitivity  Studies 

4.1  Variation  of  Damping  Parameters  of  Form  2,  Equation  2.8 

A  sensitivity  study  was  performed  on  the  model  which  provided  the  best  results 
for  roll  response,  equation  2.8,  in  which  damping  from  the  wave  excitation  was 
disr^arded  and  only  structural  damping  was  considered.  For  each  wave  test,  the 
nonlinear  damping  parameter  is  held  fixed  and  the  linear  damping  parameter  varied  by 
a  wide  range  to  determine  the  effects  on  simulated  refuse.  The  linear  parameter  is 
flien  held  fixed  and  the  nonlinear  damping  parameter  is  varied  by  a  wide  range  to 
determine  the  resulting  effects  on  simulated  roll  response. 


4.1.1  Regular  Wave,  H=  6  Ft,  T=  5  Sec  (SB26) 

The  previously  established  coefficients  for  this  case  were  a  linear  damping 
parameter  of  0.09  and  a  nonlinear  damping  parametm:  of  0.07.  Table  4.1.  shows  a 
constant  decrease  in  standard  deviation  as  damping  parametm's  increase  for  both  the 
linear  and  nonlinear  variations. 


1  Linear 

Damping 

Parameter 

Nonlinear 

Damping 

Parameter 

Standard 

Deviation 

linear 

Damping 

Parameter 

Nonlinear 

Damping 

Parameter 

Standard 

Deviation 

0.03 

0.07 

9.89 

0.09 

0.01 

9.54 

0.05 

0.07 

9.77 

0.09 

0.03 

9.48 

0.07 

0.07 

9.59 

0.09 

0.05 

9.42  1 

0.09 

0.07 

9.33 

0.09 

0.07 

9.33 

0.11 

0.07 

9.00 

0.09 

0.09 

9.28 

0.13 

0.07 

8.60 

0.09 

0.11 

9.20 

0.15 

0.07 

8.15 

0.09 

0.13 

9.12 

Table  4.1  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB26) 
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Figures  4.1a  -  4.1e  display  the  amplitude  and  phase  changes  which  result  from 
varying  the  linear  damping  parameter.  The  trend  observed  from  Figure  4.1  is  that 
with  the  lower  linear  damping  parameter,  the  re^nse  amplitude  is  greater  than  the 
measured.  In  addition,  the  phase  of  the  predicted  response  lags  the  measured  by 
approximately  0.25  secs.  Figure  4.  la.  As  the  linear  damping  parameter  proceeds 
through  the  range,  the  amplitude  deoeases  and  the  phase  shift  transitions  from  a  0.25 
sec  lag  to  leading  the  measured  reqwnse  by  approximately  0.25  secs,  Figure  4.1e. 

The  effects  of  the  variation  of  the  nonlinear  damping  parameters  is  not  as 
evident.  Figures  4.2a  -  4.2e  display  tiie  change  in  roll  response  with  variation  in  tiie 
nonlinear  damping  tmm.  The  predicted  response  phase  is  almost  identical  to  the 
measured  re^nse.  There  is  not  a  visual  change  in  the  amplitude  as  the  nonlinear 
parameters  proceed  through  the  range.  However,  a  comparison  of  standard  deviations 
shows  a  slight  decrease  as  the  nonlinear  parameter  increases,  indicating  this  term  does 
provide  a  contribution  to  the  response. 

4.1.2  Regular  Wave,  H=  6  Ft,  T=  6  Sec  (SB27) 

The  linear  and  nonlinear  damping  parameters  for  tius  case  were  estahlishfid  as 
0.10  and  0.07  respectively.  Table  4.2  indicates  a  constant  decrease  in  standard 
deviation  as  the  damping  parameters  increase.  Figure  4.3  represents  the  resulting 
changes  when  tiie  nonlinear  damping  parameters  were  held  fixed  and  the  linear 
danqiing  paramet^  were  varied.  It  is  evident  that  amplitude  decreases  with  an 
increase  in  damping  parameter.  Figure  4.4a  -  4.4e  rqiresents  the  results  of  holding 
the  linear  damping  parameter  fixed  and  varying  the  nonlinear  damping  parameter. 

The  simulated  reqxmse  amplitude  decreased  through  the  range  of  damping 
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Figure  4.1  Variation  of  Linear  Damping  Parameter,  Constant  Nonlinear  Damping 
Parameter,  =  0.07,  H=  6.6  Ft,  T=  5  Sec  (SB26,  eqn  2.8) 
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Figure  4.2  Variation  of  Nonlinear  Damping  Parameter,  Constant  Linear  Damping 
Parameter,  =  0.09,  H=  6.6  Ft,  T=  5  Sec  (SB26,  eqn  2.8) 
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paramet^s.  It  appears  the  phase  shift  remains  unchanged  throughout  the  range. 
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1  0.12 

0.07 

7.90 

0.10 

0.09 

8.43 

1  0.14 

0.07 
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1  0.16 
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0.10 
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Table  4,2  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB27) 
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Figure  4.3  Variation  of  Linear  Damping  Parameter,  Constant  Nonlinear  Damping 
Parameter,  =  0.07,  H=  6  Ft,  T=  6  Sec  (SB27,  eqn  2.8) 

--  measured,  ..  predicted 
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Figure  4.4  Variation  of  Nonlinear  Damping  Parameter,  Constant  Linear  Damping 
Parameter,  =  0.10,  H=  6  Ft,  T=  6  Sec  (SB27,  eqn  2.8) 

—  measured,  ..  predicted 
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4.1.3  Regular  Wave,  H=  7.2  Ft,  T=  8  Sec  (SB29) 

The  damping  parameters  established  for  test  case  SB29  are  0.33  and  0.40. 
There  is  a  constant  decease  in  standard  deviation  as  the  damping  parameters  increase 
as  shown  in  Table  4.3. 
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0.50 
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0.50 

0.40 

2,56 

0.33 

0.55 
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Table  4.3  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB29) 


Predicted  response  amplitude  is  greater  than  measured  when  die  nonlinear 
damping  parameter  is  held  fixed  and  the  linear  parameter  is  varied.  The  phase  of  the 
predicted  response  leads  the  measured  at  the  lower  linear  damping  parameters.  As 
the  linear  parameter  proceeds  through  the  range,  the  phase  shift  changes  and 
eventually  the  predicted  response  lags  the  measured  by  approximately  0.5  sec.  These 
results  are  diq>layed  in  Figures  4.5a  -  4.5e.  Holding  the  linear  damping  parameters 
fixed  and  varying  the  nonlinear  produces  no  significant  change  in  amplitude  or  phase 
differaice,  as  diown  in  Figures  4,6a  -  4.6e. 
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Figure  4.5  Variation  of  Linear  Damping  Parameter,  Constant  Nonlinear  Damping 
Parameter,  =  0.40,  H=  7.2  Ft,  T=  8  Sec  (SB29,  eqn  2.8), 

—  measured,  ..  predicted 
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Figure  4.6  Variation  of  Nonlinear  Damping  Parameter,  Constant  Linear  Damping 
Parameter,  =  0.33,  H=  7.2  Ft,  T=  8  Sec  (SB29,  eqn  2.8), 

—  measured,  ..  predicted 
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4.1.4  Regular  Wave,  H=  5.7  Ft,  T=  10  Sec  (SB30) 

Test  case  SB30  is  the  only  test  that  displays  sensitivity  to  change  in  damping 
parameters.  Table  4.4  shows  the  change  in  standard  deviation  over  the  range  of 
damping  parameters.  The  sensitivity  is  very  evident  in  Figures  4.7a  -  4.7e  which 
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0.01 

1.89  1 

0.02 
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0.07 

1.70 
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0.03 
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1  0.11 

0.07 

1.66 

0.03 

0.13 

1.84 

Table  4.4  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB30) 


show  the  rapid  change  in  simulated  response.  The  superharmonic  virtually  disappears 
once  the  linear  damping  parameter  reaches  0.05.  Figure  4.8  displays  the  entire  1024 
sec  time  series  for  this  case.  It  is  evident  that  the  damping  parameter  has  a  major 
effect  on  the  occurrence  and  magnitude  of  the  superharmonic.  The  standard  deviatirm 
changes  rapidly  with  the  increased  damping  parameter  until  the  siq)erharmonic  is 
eliminated  at  0.05  damping.  This  is  emphasized  by  the  fact  that  when  the  linear 
danq)ing  parameter  varies  from  0.01  to  0.05,  the  difference  in  standard  deviation  is 
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Figure  4.7  Variation  of  Linear  Damping  Parameter,  Constant  Nonlinear  Damping 
Parameter,  =  0.07,  H=  5.7  Ft,  T=  10  Sec  (SB30,  eqn  2.8) 

—  measured,  ..  predicted 


0-07  time  (seconds) 

Figure  4.8  Variation  of  linear  Damping  Parameter,  Constant  Nonlinear  Damping 
Parameter,  =  0.07,  H=  5.7  Ft,  T=  10  Sec  (SB30,  eqn  2.8) 


0.8025.  However,  when  the  damping  parameter  is  varied  from  0.05  to  0.11  the 
change  in  standard  deviation  is  only  0.08.  A  variation  in  the  nonlinear  damping 

parameter  has  very  Uttte  effect  on  the  response  ampUtude  or  superharmonic  response 
as  shown  in  Figure  4.9. 

4.1.5  Randcnn  Wave  Tests 

Tables  4.5  and  4.6  show  the  sensitivity  of  the  model  to  Bretschneider  and  white 
noise  wave  excitation.  Both  damping  models  exhibit  a  significant  change  in  standard 
deviation  over  the  range  of  linear  damping  parameters.  Li  each  case,  the  standard 
deviation  decreases  by  over  50%  as  the  linear  parameter  increases. 
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Table  4.5  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB25) 
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Figure  4.9  Variation  of  Nonlinear  Damping  Parameter,  Constant  Unear  Damping 
Parameter,  =  0.03,  H=  5.7  Ft,  T=  10  Sec  (SB30,  eqn  2.8) 

—  measured,  ..  predicted 
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Table  4.6  Variation  in  Damping  Parameters  and  Resulting  Standard  Deviations  (SB33) 


4.2  Variation  of  Nonlinear  Restoring  Moment 

All  of  die  models  considered  in  this  analysis  have  used  a  13th-order  polynomial 
restoring  moment.  A  sensitivity  test  is  conducted  to  determine  the  effects  on  the 
predicted  roll  response  whoi  lower-order  polynomials  are  used  to  fit  the  actual 
restoring  momrait  curve. 


4.2.1  Cubic  lit  to  Restoring  Moment  Curve 

Figure  4.10  shows  the  actual  restoring  moment  curve  and  the  3rd-order 
polynomial  used  to  provide  an  analytical  fit  to  the  curve.  This  cubic  term  is  then 
substituted  into  equation  2.8.  The  resulting  predicted  roll  response  is  shown  in  Figure 
4.11.  A  comparison  of  the  measured  and  predicted  spectral  densities  for  roll 
response.  Figure  4.12,  indicates  this  polynomial  is  acceptable  although  it  does  not 
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Figure  4.10  Comparison  of  Numerical  (PauUing,  1995)  and  Analytical  (Yim  et  al. 
1995)  Restoring  Momrat  Curves,  Cubic  Fit,  -  numeric^,  ~  analytical.  6135= 
3.9502e6,  B,i4=  -5.1930e6 
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Figure  4.11  Comparison  of  Measured  and  Predicted  Roll,  Cubic  Restoring  Moment 
H=  7.2  Ft,  T=  8  Sec  (SB29),  -  measured,  ..  predicted 


frequency  (Hz) 


Figure  4.12  Comparison  of  Spectral  Densities  (Deg'^/Hz),  Measured  and  Predicted 
Roll,  Cubic  Restoring  Moment,  H=  7.2  Ft,  T=  8  Sec  (SB29) 

-  measured,  —  predicted 
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provide  and  "exact"  match  of  the  magnitude  of  the  predict  roll  response  relative  to  the 
ones  with  higher-order  terms(see  following  section). 

4^^  5th,  9th  and  13th-Order  Restoring  Moment  Curves 

Figure  4.13  shows  the  5th-order  polynomial  fit  to  the  restoring  moment  curve. 
The  resulting  predicted  roll  re^nse  is  shown  in  Figure  4. 14.  A  comparison  of  die 
measured  and  simulated  spectral  densities,  Figure  4.15,  confirms  that  the  5th-order 
restoring  moment  polynomial  provides  accurate  results.  Figures  4. 16  and  4.17  show 
the  9th  and  ISth-cuda  polynomial  fits  to  the  restoring  momrat  curve  and  the  resulting 
roll  responses,  shown  in  Figures  4.18  and  4.19. 

This  s«isitivity  study  indicates  that  a  3rd  or  5th-order  polynomial  fits  die 
restoring  moment  curve  with  sufficient  accuracy  to  provide  valid  roll  re^nse 
simulation.  However,  it  should  be  noted  that  it  does  not  require  a  significant  amount 
of  time  or  effort  to  use  the  13th-order  fit,  and  this  high-order  polynomial  becomes  a 
requirement  to  accurately  predict  roll  reqxmse  at  large  roll  angles. 
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Figure  4.13  Comparison  of  Numerical  (Paulling,  1995)  and  Analytical  (Yim  et  al. 
1995)  Restoring  Momait  Curves,  5th  Order  Restoring  Moment, 

-  numerical,  -  analytical.  613^=  5.5384e6,  Bn  5=  -1.1747e7,  Bas=  5.2136e6 
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Figure  4.15  Comparison  of  Spectral  Densities  (Deg^/Hz),  Measured  and  Predicted 
Roll,  5th  Order  Restoring  Moment,  H=  7.2  Ft,  T=  8  Sec  (SB29) 

-  measured,  —  predicted 
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Figure  4.18  Comparison  of  Measured  and  Predicted  RoU,  9th  Order  Restoring 
Moment,  H=  7.2  Ft,  T=  8  Secs  (SB  29),  -  measured  roU,  ..  predicted  roU 
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5.0  Ciunparison  of  Single,  2-  and  3-  D^ree-of-FVeedom  Models 

The  final  objective  of  this  study  is  to  compare  the  results  of  roll  response  for  a 
3DOF  and  a  2DOF  model  with  the  SDOF  model.  A  comparison  of  the  5,  8  and  10 
sec  regular  wave  test  cases  is  done  as  well  as  the  Bretschneider  spectrum.  Table  5. 1 
lists  a  comparison  of  damping  parameters  for  each  of  these  cases.  It  is  evident  ftom 
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Table  5.1  Damping  Parameters  for  SDOF,  2DOF  and  3DOF  Models 


this  table  that  th^  is  a  wide  range  of  damping  parameters  among  the  different 
models.  There  does  not  appear  to  be  any  trend  established  by  the  damping 
parameters  for  each  wave  test.  In  some  instances,  the  parameters  are  higher  in  the 
SDOF  model  and  in  others  the  2DOF  and  3DOF  exhibit  higher  parameters. 

Figures  5.1a  -  5.1c,  which  show  response  to  5  sec  wave  excitation,  indicate 
that  all  three  models  simulate  the  amplitude  of  roll  response  very  well.  The  3DOF 
and  SDOF  model  match  the  phase  identically,  however,  the  2DOF  model  leads  the 
measured  response  slightly. 


Wave  ((t).  Roll  (Deg)  Wave  (II).  Roll 


(c)  SDOF  Model  time  (seconds) 


Figure  5.1  Comparison  of  SDOF,  2DOF  and  SDOF  Models,  H=  6.6  Ft,  T=  5  Sec 
(SB26) 
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All  three  models  accurately  predict  the  roll  response  for  die  8  sec  wave 
excitation  case,  Figures  5.2a  -  5.2c.  The  amplitude  as  well  as  phase  are  nearly 
idaitical  in  each  case.  It  is  interesting  to  note  that  each  of  the  models  required 
unusually  high  damping  parameters  for  this  particular  wave  case. 

Figures  5.3a  and  5.3b  compare  SDOF  and  3DOF  responses  for  the  10  sec  wave 
test  case.  This  is  the  only  case  observed  whwe  the  higher-degree-of-fineedom  model 
exhibits  better  prediction  ciqjabilities  than  the  SDOF  one.  It  should  be  noted  that  the 
3DOF  and  2DOF  models  used  analytical  input  and  the  SDOF  used  measured  input. 

The  SDOF  model  predicts  the  subharmonic  activity  seen  in  Figure  5.4b.  Figure  5.4 
compares  the  aitire  time  series  of  the  2DOF  and  SDOF  models.  Again,  the  multi- 
degree-of-freedom  model  predicts  superharmonic  response  much  more  accurately. 
However,  results  are  not  available  to  show  prediction  capabilities  using  measured 
input  for  all  three  models. 

The  final  comparison  is  for  die  Bretschneider  Spectrum  test  case.  This  is  the 
<Mily  test  case  in  which  all  models  exhibit  damping  parameters  which  are  close  in 
magnitude.  Figures  5.5a  -  5.5c  indicate  the  SDOF  and  3DOF  models  provide  an 
accurate  prediction  in  this  wave  case.  The  2DOF  model.  Figure  5.5b,  is  obsaved  to 
be  less  capable  of  producing  the  same  results. 

A  more  thorough  study  is  recommended  to  determine  the  adequacy  of  each 
model  and  to  perhaps  determine  some  correlations  betweai  the  models.  These 
preliminary  results,  however,  indicate  the  SDOF  model  may  be  equally  as  accurate 
and  in  some  cases  more  accurate  than  the  muM-d^ees-of-freedom  models.  Use  of 
the  single  degree  of  freedom  model  can  be  an  excellent  tool  for  future  research. 


(a)  3DOF  Mcxiel,  predicted  results 


(c)  SDOF  Model  time  (seconds) 

Figure  5.2  Comparison  of  3DOF,  2DOF  and  SDOF  Models,  H=  7.2  Ft,  T=  8  Sec 
(SB29) 
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(a)  3DOF  Model,  predicted  results 


(b)  SDOF  Model 


Figure  5.3  Compariswi  of  3DOF  and  SDOF  Models,  H=  5.7  Ft,  T=  10  Sec 
(SB30) 


roll  (degrees) 


85 


time  (seconds) 
(b)  SDOF  Model 


Figure  5.4  Comparison  of  2DOF  and  SDOF  Freedom  Models,  H=  5.7  Ft,  T=  10 
Sec  (SB30) 
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6.0  Concluding  Remarks 

In  this  report,  four  different  forms,  each  identical  except  for  a  variation  in  the 
damping  term  were  studied.  From  these  equations,  an  SDOF  model  to  predict  barge 
roll  motion  was  identified.  The  form  which  provided  the  most  accurate  predicted 
responses  proved  very  capable  of  simulating  barge  response  in  both  the  regular  and 
random  wave  environment  with  the  excq)tion  of  test  SB30.  A  sensitivity  study  was 
conducted  to  determined  the  effects  of  varying  damping  parameters.  In  addition,  the 
restoring  moment  term  was  examined  to  determine  the  effects  of  using  3rd,  5th,  9th, 
and  13th-order  terms.  Finally,  a  comparison  was  made  between  the  SDOF,  2DOF 
and  3DOF  models  to  identify  the  strengths  and  weaknesses  of  each  model  among  the 
different  wave  test  cases. 

6.1  Conclusions 

1)  The  model  which  provided  the  most  accurate  response  simulation  was 
equation  2.8.  This  form  assumes  the  relative  motion  effects  can  be  neglected.  As  a 
result,  only  the  barge  roll  velocity  is  considered  in  tiie  damping  terms. 

2)  The  SDOF  model  was  cq)able  of  predicting  barge  response  for  all  regular 
and  random  wave  test  cases  except  SB30,  in  which  the  10  sec  wave  excitation  was  a 
multiple  of  the  barge’s  natural  roll  frequency.  As  a  result,  the  response  motion 
contained  a  superharmonic  which  was  difficult  to  simulate. 

3)  A  constant  damping  parameter  cannot  be  used  to  predict  all  responses 
throughout  the  wave  excitation  spectrum.  The  damping  parameters  for  the  regular 
and  random  wave  excitations  summarized  in  Table  6.1,  show  no  obvious  pattern  of 
behavior. 
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Wave  Test 

Linear  Damping 
Parameter 

Nonlinear  Damping 
Parameter 

H=  6.6  Ft,  T=  5  Sec 
(SB26) 

0.09 

0.07 

H=  6  Ft,  T=  6  Sec 
(SB27) 

0.10 

0.07 

H=  7.2  Ft,  T=  8  Sec 
(SB29) 

0.33 

0.40 

1  H=  5.7  Ft,  T=  10  Sec 
(SB30) 

0.03 

0.07 

Bretschndder  Spectrum 

0.04 

0.04 

White  Noise  Spectrum 

0.07 

0.08 

Table  6.1  Summary  of  Damping  Parametm^  for  Form  2,  Equation  2.8 


4)  A  sensitivity  study  of  the  restoring  moment  term  determined  that  a  3rd  or 
5th-order  polynomial  is  sufficient  for  simulation  when  the  magnitude  of  the  roll 
response  is  not  large,  although  the  3rd-order  polynomial  tends  to  slightly  underpredict 
response  at  times.  There  is  no  discemable  differrace  in  predicted  response  provided 
by  a  9th  or  13th-order  restoring  moment  polynomial.  The  incremental  computational 
time  using  these  higher-order  polynomials  is  minimal.  In  addition,  when  roll 
re^nse  becomes  large,  the  higher-order  polynomials  become  a  requirement. 

5)  When  comparing  the  SDOF  model  with  the  multi-degree-of-freedom 
models,  the  results  are  mixed.  In  the  regular  wave  tests,  the  SDOF  model  provides 
the  same  results  as  the  3DOF,  except  in  the  10  second  wave  test  case.  The  ability  of 
SDOF  and  3DOF  models  to  simulate  random  response  is  remarkably  better  th^n  the 
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2D0F  model.  These  observations  combined  with  the  significant  time  saving  makes 
the  SDOF  very  desirable  to  use. 

6.2  Recommendations  for  Future  Studies 

1)  Further  studies  in  variation  of  the  damping  terms  may  improve  the  ability 
to  predict  superharmonic  response.  Varying  the  effects  of  relative  motion,  roU 
velodty  and  wave  dope  velocity  could  possibly  produce  better  results. 

2)  A  computer  program  is  currently  under  development  to  predict  the 
probability  of  barge  cs^size.  A  comparison  of  the  SDOF,  2DOF  and  SDOF  models 
can  be  used  in  this  program  to  establish  further  strengths  and  weaknesses  of  each 
model.  Although  simplistic  compared  to  the  2DOF  and  SDOF  models,  the  SDOF 
model  may  provide  equally  adequate  results  with  much  less  time  invested. 

S)  In  determining  the  damping  parameters  in  this  study,  the  standard  deviation 
of  the  amplitude  was  used  as  the  criterion  for  determination.  This  assumes  the  phase 
shift  is  unimpOTtant.  The  focus  of  an  alternate  study  could  be  to  use  a  least-square 
criterion  which  weighs  bofli  amplitude  and  phase  shift  in  the  error  calculation. 

4)  In  the  random  wave  cases,  a  comparison  of  spectral  densities  showed  a 
shift  in  peak  frequencies  between  the  measured  and  predicted  spectrum.  Variations  in 
input  data  such  as  dominant  period  of  wave  spectrum  could  be  analyzed  to  determine 
if  this  shift  may  be  eliminated  or  decreased. 


References 


90 


Bartel,  W,  Modelling,  Validation,  and  Simulation  of  Multi-Degree  of  Freedom 
Nonlinear  Stochastic  Barge  Motions,  1996 

Chakrabarti,  S.K.,  Hydrodynamics  of  Offshore  Structures,  1987. 

Falzarano,  J.,  Shaw  S.  W.  and  Troesch,  A.  W.,  Application  of  Global  Methods  for 
Analyzing  Dynamical  Systems  to  Ship  Rolling  Motion  and  Capsizing, 
ImemationalJoumal  ofBijurcation  and  Chaos,  1992,  2,  101-115. 

McAllister,  M.,  "Evolution  and  State  of  the  Art  of  the  Log  Barge,"  Marine 
Technology,  SNAME  1995,  32,  No.2  ,  132-139. 

Paulling,  J.R.,  "Capsize  Assessment  Technology  Program,  Comparison  of  Model 
Tests  and  Results  of  Numerical  Simulations,"  Report  Submitted  to  Naval 
Facilities  Engineering  Service  Center,  Port  Hueneme,  CA,  Jan  1995. 

Soliman  M.  and  Thompson,  J.M.  T.,  Transition  and  Steady  State  Analysis  of  r.apsi7ft 
Phenomena,  Applied  Ocean  Research,  1991, 13,  82-92. 


Yim,  S.C.S.,  Bartel,  W.A.,  Huang,  E.T.,  Lin,  H.,  "Nonlinear  Roll  Motion  and 
Capsizing  of  Barges  in  Random  Seas,  Part  I:  Modelling  and  System 
Identification",  Report  Submitted  to  Naval  Facilities  Engineering  Service 
Center,  Port  Hueneme,  CA,  1995,  1-41. 


